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ABSTRACT 

We investigate the behavior of the fifth force in voids in chameleon models using the 
spherical collapse method. Contrary to Newtonian gravity, we find the fifth force is 
repulsive in voids. The strength of the fifth force depends on the density inside and 
outside the void region as well as its radius. It can be many times larger than the 
Newtonian force and their ratio is in principle unbound. This is very different from 
the case in halos, where the fifth force is no more than 1/3 of gravity. The evolution 
of voids is governed by the Newtonian gravity, the effective dark energy force and 
the fifth force. While the first two forces are common in both ACDM and chameleon 
universes, the fifth force is unique to the latter. Driven by the outward-pointing fifth 
force, individual voids in chameleon models expand faster and grow larger than in a 
ACDM universe. The expansion velocity of the void shell can be 20% to 30% larger 
for voids of a few Mpc/h in radius, while their sizes can be larger by ~ 10%. This 
difference is smaller for larger voids of the same density. We compare void statistics 
using excursion set theory; for voids of the same size, their number density is found 
to be larger in chameleon models. The fractional difference increases with void size 
due to the steepening of the void distribution function. The chance of having voids 
of radius ~ 25 Mpc/h can be 2.5 times larger. This difference is about 10 times 
larger than that in the halo mass function. We find strong environmental dependence 
of void properties and population in chameleon models. The differences in size and 
expansion velocity with GR are both larger for small voids in high density regions. 
In general, the difference between chameleon models and ACDM in void properties 
(size, expansion velocity and distribution function) are larger than the corresponding 
quantities for halos. This suggests that voids might be better candidates than halos 
for testing gravity. 
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1 INTRODUCTION 

Models of modified gravity (MG) are introduced to explain 
the observed accelerating cosmic expansion, without invok- 
ing a cosmological constant in the Einstein equation. Scalar- 
tensor gravity theories are among those that are well re- 
ceived recently. In these theories, the scalar field is cou- 
pled to matter, triggering an extra fifth force which leads 
to an universal enhancement of gravity. The enhanced grav- 
ity violates existing robust tests of general relativity (GR) 
in the solar system, so that only theories with a screening 
mechanism to suppress the fifth force in high density regions 
are observationally viable (e.g. Khoury fc Weltman 2004| ). 
Gravity is therefore back to GR in the early universe, as well 
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as in the vicinity of virilized objects where the local density 
is sufficiently high. MG models like chameleon gravity can 
therefore pass the tests of current constraints from the so- 
lar system ( Khoury fc Weltman||20d4 |. Nevertheless, struc- 
ture formation in these models should be somewhat differ- 
ent from that of the standard A-cold-dark-matter (ACDM, 
where A represents the cosmological constant) paradigm. In 
low density regions of the universe, the fifth force is weakly 
or not suppressed, so that dark matter and ordinary mat- 
ter will feel this extra force and hence evolve differently from 
the GR case. Qualitatively, one may expect structure to form 
earlier in MG than in GR with the help of enhanced grav- 
ity. Indeed, halos are found to be more massive and more 
abundant in simulations of f{R) gravity ( Li et al^|2012 \ as 
compared to GR at the same epoch. Similarly, voids appear 
to be larger and emptier in MG. These qualitative results 
seem to point in the same direction as some recent observa- 
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tional facts, which have been shown to be in tension with a 
ACDM universe. 

Firstly, some galaxy clusters detected using X-ray and 
lensing techniques at high redshift are found to be too mas- 



sive and have formed too early (e.g. Enqvist et al. 2011 
Jee et al.|20Tl| |Hoyle et al.|2011[|Holz fc Perlmutter|2012 1 

The probability of the existence of those massive clusters 
in ACDM is prohibitively small, but see Hotchkiss (2011!); 



Harrison fc Coles| p()12| ); [Hoyle et al.| ( |2012[ ); [Waizmann et| 
al.| ( |2012a|b| |. Introducing non-Gaussianity can ease this ten- 
sion, but the /nl parameter required to fit the data is usu- 
ally too high, which is in tension with other observational 
constrains like the cosmic microwave background (CMB). 



Secondly, the detected integrated Sachs Wolfe (ISW) ( Sachs 
fc Wolfe|l967 1 signal from the stacking of 4-deg^-size regions 



of the CMB corresponding to the SDSS super clusters and 
super voids is found to be 2 — Scr higher than estimations 



from simulations (Granett et al. 2008 Papai et al. 20111 



This tension with the ACDM paradigm is perhaps more than 
3 — (T as suggested in Nadathur et al. ( 2012 1. Accounting for 



non-linear ISW effect by using simulations of full-sky ISW 
maps from ( Cai et al.|2010 1 , the tension remains nearly un- 
changed ( Flender et al.|2012 l. Similar conclusions are found 



by independent study of ( Hernandez-Monteagudo fc Smith 



2012 1. If one assumes that the expansion history of the uni- 
verse is given by the concordance ACDM model, then one 
plausible explanation of this discrepancy is that the abun- 
dance of structure in the real Universe may be greater than 
expected, i.e., there might be more clusters and super clus- 
ters, and voids might have grown larger and deeper. This ex- 
planation seems to coincide with the first tension mentioned 
above. Again, one could perhaps use this data to constrain 
non-Gaussianity, and find a large /nl, but an alternative 
solution might be to modify gravity. 

In this work, we explore the difference of structure for- 
mation in GR and chameleon models of MG. Using the 
spherical collapse model and excursion set theory [Bond et aL] 
(19911, we investigate individual void properties and the 
void volume distribution function in these two models. We 
also compare the relative merits of distinguishing between 
GR and MG using voids or halos; predictions for the latter 



have been addressed by Li & Efstathiou (20121 



One common way to distinguish MG from GR is by 
looking at the difference between the lensing mass and dy- 



namical mass of halos ( 


Feix et al.|2008, ,Schmidt.2010| Jain 


fc Khoury 


20101 |Zhao et al.|2010||Clampitt et al.|2012||Lam 


et al.|2012 


\. The chameleon model studied here predicts that 



such a difference is at most 1/3 between the screened and 
unscreened cases. At present, it is still very difficult to have 
mass estimates of halos which achieve this level of accu- 
racy, partly due to the difficulty of looking for unscreened 
objects. To realize the 1/3 difference, such objects must be 
both small, so that they are not self-screened, and located 
in low density environments, so as not to be screened by the 
environment. Voids, however, are usually very low in density 
so that the fifth force is unscreened inside. Furthermore, we 
show that the strength of the fifth force may be relatively 
stronger than that of Newtonian gravity in voids. This may 
lead to a larger difference of void properties from GR than 
that of halos. 

The outline of this paper is as following: In section [2j 
we give a brief summary of the coupled scalar field gravity, 



of which the chameleon model is an example. In section [3j 
we solve the scalar field profile for voids in this model and 
highlight interesting differences of the fifth force to Newto- 
nian gravity in voids. In section [4j we extend the spherical 
collapse model to solve for the evolution of shells in voids in 
this model and identify the best regimes to distinguish this 
model from GR. Section [5] presents the first crossing barrier 
for voids, and incorporates the moving barrier and environ- 
mental dependence of void formation to the excursion set 
theory to calculate a void volume distribution function. We 
summarize our results and consider possible ways to test 
MG in voids in section [G] 



2 THE CHAMELEON THEORY 

This section lays down the theoretical framework for investi- 
gating the effects of a coupled scalar field in cosmology. We 
shall present the relevant general field equations in 
and then specify the models analyzed in this paper in 



2.1 
T2I 



2.1 Cosmology with a Coupled Scalar Field 

The equations presented in this subsection can be found inO] 



fc Zhao^ ( .2009, ,2010j ) ; |Li fc Barrow] \2Q11\ , and are presented 
here only to make this work self-contained. 
We start from a Lagrangian density 

£ = ^ [Ml^R - V>Va0] + V{<l)) - C(0)(£dm + £s), (1) 
in which R is the Ricci scalar; the reduced Planck mass is 



Mpi = I/VSttG with G being the gravitational constant; 
and £dm and Lq are respectively the Lagrangian densities 
for dark matter and standard model fields. is the scalar 
field and V((f>) its potential; the coupling function C{(j)) char- 
acterises the coupling between and matter. Given the func- 
tional forms for V{(f)) and C{(f)), a coupled scalar field model 
is then fully specified. 

Varying the total action with respect to the metric gab, 
we obtain the following expression for the total energy mo- 
mentum tensor in this model: 



Tab = ^a(l)^b4> - dab 



\/''\/c(l)-V{(t>) 



+ cmTT + T!b), 



where T^b^ and Tab are the energy momentum tensors for 
(uncoupled) dark matter and standard model fields. The ex- 
istence of the scalar field and its coupling change the form of 
the energy momentum tensor, leading to potential changes 
in the background cosmology and structure formation. 

The coupling to a scalar field produces a direct in- 
teraction (fifth force) between matter particles due to the 
exchange of scalar quanta. This is best illustrated by the 
geodesic equation for dark matter particles 

where r is the position vector, t the (physical) time, <I> the 
Newtonian potential and V is the spatial derivative; = 
dC/d</i. The second term on the right hand side is the fifth 
force, with potential lnC(</)). 

To solve the above two equations we need to know both 
the time evolution and the spatial distribution of 0, i.e. we 
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need the solutions to the scalar field equation of motion 
(EOM) 



(4) 



where p = pDM + pb, the sum of dark and baryonic matter 
densities. Equivalently 



where we have defined 



(5) 



(6) 



The background evolution of can be solved easily given the 
present-day value of p since p cc a~^. We can then divide (j> 
into two parts, <j> — (j> + 5(f), where (j> is the background value 
and 54> is its (not necessarily small nor linear) perturbation, 
and subtract the background part of the scalar field equation 
of motion from the full equation to obtain the equation of 
motion for 5<j). In the quasi-static limit in which we can 
neglect time derivatives of 5(f) as compared with its spatial 
derivatives (which turns out to be a good approximation on 
galactic and cluster scales), we find 



dC(<^) 
A(t) ' 



dC(,^)_^dV(0) dV(0) 



d(t> 



d(f) 



(7) 



where p is the background matter density. 

The computation of the scalar field (fi using the above 
equation then completes the computation of the source term 
for the Poisson equation 



[ptot + 3ptot] 
[pCi(^)-2Vm, 



(8) 



where we have neglected the kinetic energy of the scalar field 
because it is always very small for the model studied here. 



force is strongly suppressed in high density regions where 
(f> acquires a large mass, 2> (H is the Hubble ex- 
pansion rate), and thus the fifth force cannot propagate far. 
The suppression of the fifth force is even stronger at early 
times, and thus its influence on structure formation occurs 
mainly at late times. The environment-dependent behaviour 
of the scalar field was first investigated by |Khoury fc Welt-| 



( 2004 1, and is often referred to as the 'chameleon effect'. 



3 STATIC UNDERDENSITY SOLUTIONS 

The radial profile of a chameleon-type scalar field has been 
studied in detail for spherical overdensities, in which cases 
a simple analytical formula for the fifth force has been de- 



rived ( Khoury & Weltman 2004 1 and shown to agree well 
with the numerical simulations ( |Li et aL||2012[ ). We know 
from these previous studies that, depending on its size and 
environment, a spherical overdensity could develop a thin 
shell which is a region of fast change of (f){r) with respect to 
r, and approximately only the matter contained in this shell 
contributes to the fifth force on a particle at the edge of the 
overdensity. If the shell is thin the fifth force is much weaker 
than gravity (the latter coming from all mass contained in 
the overdensity), while if its thickness becomes comparable 
to the radius of the overdensity, the fifth force approaches a 
constant ratio to gravity. For our fiducial model this ratio is 
27^ and we choose the coupling 7 such that 27^^ = 1/3, so 
that the maximum deviations from GR match those of f{R) 
models. 

Unfortunately, no analytical approximation for the fifth 
force is known for the case of underdensities. It is our task 
in this section to study (fi{r) in underdensities and the fifth 
force which results. We will see that the maximum ratio of 
27^ = 1/3 will no longer apply in this case: in voids the fifth 
force can have much stronger effects than gravity. 



2.2 Specification of Model 

As mentioned above, to fully fix a model we need to specify 
the functional forms of V{(f>) and C{(j)). Here we will use the 
models investigated by |Li fc Zhao| ( |2009[ |2010[ ); |Li] ( |2011[ ), 
with 



(7(0) = exp(7<j!>/Mpi), 



and 

V{(f))^ 



PA 



e^pi-(f)/Mpi)Y 



(9) 



(10) 



In the above pA is a parameter of mass dimension four and 
is of order the present dark energy density {(f> plays the role 
of dark energy in this model) . 7, a are dimensionless pa- 
rameters controlling the strength of the coupling and the 
steepness of the potential respectively. 



We choose a < 1 and 7 > as in Li fc Zhao (2009 



20101, which ensure that 14ff((/)) has a global minimum close 
to = and that d^VEa((fi)/d(f>'^ = m| at this minimum 
is very large in high density regions. There are two con- 
sequences of these choices of model parameters: (1) (f> is 
trapped close to zero throughout cosmic history so that 
V{(f>) ~ Pa behaves as a cosmological constant; (2) the fifth 



3.1 Voids in Newtonian Gravity 

Consider a spherically-symmetric underdensity defined by 
radius r and inner and outer densities, pin and pout, such 
that Pin < Pout . First we review the forces around such voids 
in Newtonian gravity. Since C{(j>) « 1, the first term on the 
right-hand side of Eq. (|8| can be integrated once to give the 
force per unit test mass 



where 



_ GM{< x) 
Fn{x) —2 



M{<x)=^n rdx'x'poix')- 
Jo 



(11) 



(12) 



We are interested in the simplest model of a void, with top- 
hat density profile 



Poix) = 



Pin 
Pout 



for 
for 



X>r 



(13) 



(We use the notation r for the void radius and x for the ra- 
dial coordinate for the sake of continuity with later sections 
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of the paper.) The resulting force on the mass sheU at r is 
AttG 



3 



Pinr 



(14) 
(15) 



Only mass within the radius r contributes to the force on it - 
test masses inside completely empty voids where pin = feel 
no force since the pull from all the mass elements outside the 
void cancel perfectly. This is a standard, although counter- 
intuitive, result of Newtonian gravity. If pin is nonzero, the 
force on the shell is equal to that of a point particle of mass 
M(< x) which is located at x = 0, and the force is attrac- 
tive. 

Similarly, since V{(j)) ~ pA, the second term on the 
right-hand side of Eq. (|8| gives the effective force due to 
the scalar field potential (or equivalently, the cosmological 
constant). 



PAT 



(16) 



This contributes an effective repulsive force at late cosmo- 
logical times, which we call the dark energy force in this 
paper. 

3.2 Voids in Chameleon Theories 

The total force on a test particle is the sum of the New- 
tonian force, effective force from the dark energy, and the 
scalar-mediated fifth force. We will see that the fifth force 
is always repulsive in voids, in the sense that the force on 
a test mass pushes it away from the center of the void, to- 
wards the nearest wall; it aids the dark energy in emptying 
the void of matter. 

The second term on the right-hand side of Eq. ( 3| gives 
the fifth force, which for our choice of C{(j>) in Eq. (91 is 

We define the ratio of fifth to Newtonian forces as 
F5 _ 67A/P1 d(j> 
Fk rpin dx ' 



V 



(17) 



(18) 



which is constrained to be ^ 27 for overdensities. Thus 
the problem of finding the force deviations on a test particle 
in the void has been reduced to obtaining the scalar field 
profile (/I>(x)- Before solving Eq. ([5| to obtain the profile, we 
note some properties of this scalar field model which will 
simplify the solution. 

At fixed density po, our theory has an effective potential 

^cff(0) = 7:^ + PO exp(7<,!./Afpi). (19) 

[1 - exp (-0/i\/pi)J 

Call (po the field value which minimizes this potential for the 
given density. Using the facts that a <^ 1 and 0o/Mpi ^ 1 
(Sec. 2.2 1, we set dVeff/d<j) = and expand in the small 



parameter cj!>o/Afpi to find 

(j>o/Mpi 



a PA 
7 po 



(20) 



If mo is the mass of small fiuctuations about this minimum, 
then 



(21) 



so that the associated Compton wavelength Ao 
the field value at the minimum are related by 



and 



(22) 



The above analytic relations between the density and asso- 
ciated field value and compton wavelength, namely 1/po oc 
00 oc Ao, are not a general feature of scalar-tensor theories 
of gravity, nor even of chameleon models. For example, in 
the f{R) model of|Hu fc Sawicki] (|2007[l, the relation be- 



tween these three quantities has no closed form solution. 
While these analytic relations are useful in themselves, we 
now show how they can be used to decrease the void param- 
eter space from three to two variables, while simultaneously 
removing dependence on the theory parameters a and 7. 

Naively, any top-hat void of radius r and density pin 
in a uniform background density pout is dependent on three 
length scales: r, Ain, and Aout. However we show that since 
(p/Mpi <^ 1, the Planck scale drops out of the equation of 
motion, giving us the freedom to rescale the solution by one 
of these lengths. This reduces the problem to two non-trivial 
degrees of freedom. The equation of motion Eq. ([5| is given 
by 



d> _^ 2 d<j!> 



(23) 

where po(x) is again the top-hat profile of Eq. ( 13 1. Expand- 
ing to lowest order in (p/Mpi and using a 1 we have 



dV _^ 2 d<^ 
dx^ X dx 



apA 



4>o(x) <t> 



apA 

pout 



(pout \(po(x) 



(24) 



(25) 



Defining the dimensionless field ip 
yields 

d^V ^ 2 d^ _ 1 



and using ( 22 1 



^ 2d^ 

dx^ X dx 



(26) 



Then defining a dimensionless radial coordinate r = x/^out, 
the equation further simplifies to 

dV , 2d-i/) (Pout 1 



^ 2d^ 

dr^ T dr 



Mr) ^ ■ ^^^^ 
Now, from the three length scales we can form two ra- 
tios r/Aout and Aout/Ain (note that for voids we must have 
5S Aout/Ain < 1) and recast the EOM in terms of these. 
The problem is then reduced to solution of the differential 
equation 



dV 2 dip 
dr^ T dr 



1 



Aout/ Ain 
1 



for 



for 



r < r/Aout 
T > r/ Aout 



with boundary conditions 
= 0, 



dip 
d7 



iP{t — >■ 00) — 1 . 



(28) 



(29) 



Rewriting Eq. ( 18 1 in terms of the new variables and using 
Pa ~ ^APc we find 



ry(r) = 67 



Q^A A/piy^p^ 



rpin 



dip 

d7 



(30) 



-=r/Ao 



where pm is the background matter density today. 

Before describing the resulting solutions of the scalar 
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field and fifth force for reafistic underdensities, we make 



some comments about the relevance of Eq. ( 28 1 for our re- 



sults in Sections |4] and [5] Since there is no known analytical 
approximation for in underdensities, as there is in the over- 
dense case, it will be necessary to solve numerically the EOM 
at each time step for an expanding void in Sec. |4] Further- 
more, in order to obtain the void-formation barriers of Sec.[5j 
we must calculate the trajectories of many such expanding 
voids of different initial sizes and densities. While a top-hat 
underdensity intrinsically has three degrees of freedom, pin, 
Pout, and r, we have shown that two ratios formed from these 
quantities are sufficient to solve the EOM. Thus, the most 
difficult numerical challenge of Sections[4]and[5]can be over- 
come with a single two-dimensional table of dtp/dr values, 
where the derivative is evaluated at the border of the void. 



Furthermore, our recasting of Eq. ( 23 1 as Eq. ( 28 \ has no ex- 



plicite dependence on the theory parameters a and 7. Thus, 
this same 2-D table serves to calculate the void-formation 
barriers under variations in a and 7, as in Sec. |5.6| 

3.3 Radial Profile of the Scalar Field 

Now we consider the results for the radial profile of the scalar 
field in various underdensities, paying special attention to 
the value of the derivative at the void border. The left panel 
of Fig. [1] shows the dependence on void radius, r. If r is 
small, then the underdensity can be considered as a small 
perturbation on the environment and the scalar field value 
inside is very close to its value at the boundary. As r in- 
creases, however, there is increasing space for (j) to evolve 
away from the exterior value (here (pout ~ 0.8 x 10~^Mpi) 
as x/^ decreases, and therefore the scalar comes closer to 
reaching the value which minimizes the interior effective po- 
tential. Since (?!)in(out) oc l/pin(out) and Pout/pin = 5 in the 
figure, we know (jjin = 50out and see that even 160Mpc//i is 
not enough space for the scalar field to attain its minimum 
at the center of the void. 

Fig. [1] shows the dependence of the scalar field profile 
(j){r) on interior density pin in the central panel, assuming 
an exterior density equal to the cosmic mean. Here the field 
does not experience much change between the outside and 
inside of the void, growing by only 25% in the most extreme 
case, Pin = 0.1 pm- As a result, the derivative of the scalar, 
and therefore the fifth force, at the void border x/*" = 1 
must be small. However, we will see that in order to get a 
full picture of the forces involved it is necessary to consider 
the gravitational force and dark energy force as well. For 
this void the magnitude of the fifth force is about twice as 
large as Newtonian gravity, so that even this slowly varying 
4> profile results in a force that is stronger than _Fn. 

Finally, the dependence on pout is shown in the right 
panel of Fig. [l] The variations here appear more drastic, 
since only in this panel is the limiting value ^out changed 
from one curve to another. With fixed interior density, a 
denser environment for the void results in a larger change 
in the scalar and correspondingly higher derivative d<^/dx. 
Note also that due to Birkhoff's theorem, changes in pout 
do not affect the gravitational force inside the void, nor is 
the dark energy force is affected. So only from this panel 
can we infer directly that larger gradients of (j> imply greater 
deviations from GR. 

There are some interesting differences from the over- 



dense case. Consider an overdensity and underdensity each 
embedded in the same environmental density p'^ut , with cor- 
responding minimum (jiout- For the overdensity, we know (j) 
decreases from ^^ut as we move towards the center; how- 
ever by the shape of the effective potential Eq. (191, (j> is 
strictly positive, so < (l>'h-i < 4''out- The maximum change 
is therefore A(j) — </>out, no matter how great is the inte- 
rior density. In contrast, for the underdensity, (j> increases 
from 0out as we move towards the center so that has 
no such bound: pi„ can be infinitely small in principle. For 
concreteness consider the lowest curve on the right panel 
of Fig. [1] here (/!>out ~ 10~^Mpi so that for an overdensity 
A(f> < 10~®Mpi, while for the pictured underdensity we see 
Since the fifth force is proportional to 
the derivative of (j) at the void border, we expect this lack of 
upper bound on A0 for underdensities to show itself in the 
force. We turn our attention next to these results. 



3.4 The Fifth Force 

In the top panels of Fig. [2] we show the force deviation 
77 — F^/Ft^ with variations in the three physical parame- 
ters which define a void, r, pin, and pout. The first interest- 
ing feature is that rj is always negative. The fifth force in 
voids is repulsive, always pointing at the opposite direction 
of normal gravity. This is the direct consequence of the scalar 
field profile we have shown in Fig. [T] whose slope is always 
negative at the edge of the underdensity. Intuitively, this re- 
pulsion occurs due to the Yukawa potential {e~^'''^ /x) oi the 
scalar: at distances of the order of the Compton wavelength, 
the potential falls off more strongly than 1/x. Mass elements 
on the far wall of a large void are unable to cancel the pull 
of the near wall. Furthermore, even if pin is nonzero, the in- 
tegrated mass inside the shell is unable to compete with the 
denser nearby wall, and the force is again repulsive. 

Secondly, as we anticipated in Sec. |3.3| the unbounded- 
ness of the field, along with the result of Eq. ( |15[ ) that Fn 
vanishes as pin — )■ or r — >■ 0, leads to deviations which do 
not share the bound oi \ ri\ ^ 27^. Thus the relative strength 
of the fifth force can be much larger than Newtonian gravity, 
as seen in all top panels of Fig. [2] Even for common voids 
with ratio of densities pin/pout = 0.2 we can have \ri\ « 1/3, 
already reaching the upper bound for overdensities. If the 
ratio decreases to the percent level, then rj ~ —2 for the 
smallest voids. 

The left panel of Figs.[l]and top-left panel of Fig.[2]both 
show variations with respect to void radius r. Comparing 
these, we see that while the change in (p and therefore the 
fifth force increases with void radius, the deviation 77 gets 
smaller. Thus we infer that Fn increases more quickly than 
F5 in these cases due to the increasing mass enclosed within 
the larger void radius. 

In contrast, comparing the middle panel of Figs, [l] and 
top-middle panel of |2] shows that under variations in pin the 
changes in the fifth force dominate the dependence of rj. The 
net effect of decreasing the interior density is to strengthen 
the fifth force relative to gravity. 

The variations of pout in the top-right panel of Fig. [2] 
leave _Fn unaffected, so here changes in straightforwardly 
reflect changes in F5. We can unify the results of varying 
Pout and Pin by noting that increasing the density contrast 
Pout/pin generally increases the deviation from GR. 




Figure 1. Left panel: Radial profile of the scalar field in a spherical top-hat underdensity for fixed values of pj^ = 0.2 p^, pout = Pm 
and different radii r. Center panel: The same, but for fixed values of pout = pm, r = 20Mpc//i and different inner densities pin. Right 
panel: The same, but for fixed values of r = 20Mpc//i, = 0.2 pm and different outer densities Pout- Note that the horizontal axis is 
scaled with respect to void radius r, so x/** = 1 is the edge of the spherical underdensity; also we evaluate the cosmic mean density at 
the present day, Pm(^) = Pm(0). 



In principle this unboundedness of the force ratio r] in 
underdensities looks very promising for distinguishing be- 
tween GR and chameleon models. However, at late times 
when pm and pA are comparable, the repulsive dark energy 
force can dominate over Newtonian gravity where the den- 
sity is low. Fa is common in both GR and MG models but 
negligible for halos where the local density is much greater 
than the cosmic mean. The evolution of voids in MG models 
are therefore affected by F^, Fa and Fn- 

Bottom panels of Fig.[2]show the fractional difference of 
total force between MG and GR, F^/^Fn+Fa)- Comparing 
them with the top panels, we find the following: A.) like ry, 
the fractional difference decreases with radius (bottom-left) 
and increases with pout (bottom-right). This is because the 
additional Fa term is just a constant at a certain epoch. B.) 
Fc,/ [Fm+Fa) can be positive or negative, depending on the 
relative amplitude of _Fjv and Fa. The transition occurs at 
Pin = 2pA when Fn is canceled out by Fa, and the evolution 
of the system is only governed by F^. Note that the sign 
switch in Fi,/{Fn+Fa) is no more than an indicator for the 
switch of the relative strength between Fn and Fa. The 
forces Fa and F5 are always repulsive, and act to accelerate 
the expansion of void. C.) When pin is close to 2pA, the 
fractional difference can be very large. 

In summary, if we track the evolution of a spherical 
underdensity with the radius of r, in the early universe it 
is dominated by Fn, the amplitude of which decreases with 
pin - Later, the repulsive dark energy force Fa from the back- 
ground scalar field emerges to cancel part of Fat , and it helps 
to accelerate the expansion of void shells. In the mean time, 
F5 appears from the coupling of the scalar field with mass, 
and is also repulsive in voids. As the void keeps emptying 
itself, F5 becomes larger and Fa also grows with time as 
f^A increases. The amplitude of the positive Fjv+Fa keeps 
decreasing until pin = 2pA, then Fn + Fa switches sign and 
the amplitude starts increasing. F5 should also keep increas- 
ing with time as pin decreases faster than its environment 
density, which makes the density contrast inside and and 



outside the void grow larger. Overall, F5 should help to ac- 
celerate the expansion of void. In the next section, we will 
quantify this effect. 



4 EVOLVING INDIVIDUAL VOID 

With the solution of the fifth force in underdense regions, we 
can apply it to solve the equations that govern the evolution 
of a spherical underdensity in a given environment specified 
by its density. We will explore how the evolution of voids or 
underdense regions are affected by the fifth force. 



4.1 Evolution of Environment 

We have shown in the previous section that the profile of 
the scalar field and hence the fifth force depends on the lo- 
cal density as well as the density of its environment. This is 
one distinct feature of chameleon models. We therefore need 
to follow the evolution of the environment properly in order 
to calculate the fifth force. The environmental dependence 
in chameleon models has been discussed by |Li fc Efstathiou] 
(20121) and iLi & Laml (120121) for halos. We shall adopt the 



same idea of taking the environment as a spherical region 
with radius much larger than the underdensity in considera- 
tion. The exact choice of the environment size will be spec- 



ified where it is used later for the void statistics (Sec. 5.31. 
Note that for the purposes of single-shell evolution which 
we describe in this section, the environment is completely 
specified by its density relative to the cosmic mean. 

To track the non-linear evolution of the environment, 
we denote its physical radius at time t by r{t), its initial co- 
moving radius by R, and define q{t) = a(t)R. The evolution 
equation for r{t) is 



(p-2pa), 



(31) 



where p = 3M/47rr^ is the matter density in the spherical 




Figure 2. Top-left panel: Variations of the force deviation rj with underdensity radius, r. The exterior density is fixed to the cosmic 
mean today, pout = Pm- Various values of interior density are shown, with pj^ decreasing from top to bottom. Top-center panel: 
The same, but for continuous variations of pjn, fixed pout = Pm and various values of radius r, with r decreasing from top to bottom. 
Top-right panel: The same, but for continuous variations of pout, fixed pi^ = 0.1 pm, and various values of radius r, with r decreasing 
from top to bottom. Note that in all the panels, we evaluate the cosmic mean density at redshift one, pm{z) = p,7i(l). Bottom-left, 
-middle and -right panels are the same as the top-left, -middle and -right panels, but showing the fractional difference of the total force 
between MG and GR theories, Fs/I^Fn+Fa)- 



region of the environment and the constant pA ~ V{<j)) is the 
effective dark energy density. Note that Eq. ( |31[ ) assumes 
that the environment is unaffected by the fifth force. We 
make this approximation since the environments are very 
large in size and therefore the effects of the fifth force on 
them are minimal. Let us define y{t) = r{t)/q{t) and change 
the time variable to N = ln(a); derivatives with respect to A'^ 



are denoted by y' — Ay/ AN. By using Eq. (31 i, q{t) oc a{t) 



and the Friedman equation H / Hq = f2„ 



JIa, we find 



2--Q.m{N) 



y + 



{y-'-l)y = 0, (32) 



a non-linear 

'^+nA) 



equation, where Qm{N) 
and nA{N) = nA/i^me~^^ 



+ 



which 

At very early times we must have y ~ 1 and so can 
write y = 1 -\- 1 with |e| <IC 1. Substituting this into Eq. (321 
to get the linearised evolution equation for e, we find that 
e oc D+, in which D+ is the linear growth factor governed 



by the equation 



2--^m{N) 



D'+ - -Qm{N)D^ 



0, 



(33) 



and the proportionality coefficient can be found using mass 
conservation: y^{l -\- 5i) = 1 ^ e = —5i/3 oc -D+ (here Si 
is the linear density perturbation at the initial time). As a 
result, the initial conditions for y are y{ai) = 1 — 5i/3 and 
y'{ai) = -5j/3- 



Eqs. (32 33 1 , associated with their corresponding initial 



conditions, completely determine the necessary dynamics in 
the ACDM model used for the environment shell. In what 
follows we shall use j/cnv to denote the y for the environment, 
in contrast to that for the underdensity, which we shall de- 
note by yv We will reserve r for the physical radius of the 
underdensity, matching the notation of Sec. [3] 



8 J. Clampitt et al. 



4.2 Evolution of Underdensity 

The only difference between the evolution of an underdensity 
and that of its environment is the effect of the fifth force. To 
calculate the fifth force at each time-step we use a spherical 
top-hat profile, 

Pv for x^r 
^Ponv for X> r 

We assume there is no shell crossing, so that to study the 
evolution we only need to understand the motion of the shell 
at the edge. Note that this is not strictly true: for a model 



Pix) 



different than ours, Martino & Sheth ( 2009 1 have shown that 



modified gravity can cause an initially top-hat underdensity 
to have a slight density gradient near the edge. We find a 
similar effect, but it is quite small and it is beyond the scope 
of this paper to self-consistently track the deviations of the 
density profile from the top-hat. 

Denoting the density inside the underdensity by pv and 
using mass conservation, we can show that 

3 /- -3\ / r)\3 



Pv = 



= IPmtS 



1 i<^Rr 



Pm (aj/v) 



(35) 



where pm is the background matter density today. Simi- 
larly, the matter density in the environment, Ponv, can be 
expressed in terms of j/cnv as 

-3 



Pm {aye 



(36) 



Using these relations we can rewrite Eqs. ( 15 I and (17 1 in 
terms of the variables j/v and ycnv, yielding 
1 



where 



6M, 



2-Pm (ayv) 



■R 



in 



{HoR) (ay,) 
d {(t>/Mj.i) I 

dx \x^r 

dr It: 



' Ho, 



1 



=>-/Ao„ 



V Ml A 



HoR, 



(37) 



(38) 



(39) 



and tp and r are defined as in Sec. [3] The fifth-force-to- 
gravity ratio is then 



V = 



, (HoR) (ay.)-'' 



(40) 



and the evolution equation of the underdensity becomes 
r 1 



[p.{l+Tj)-2pA] 



Rewriting using y^ we obtain 



J/v + 



Vv + 



(41) 



[y-^(l+7?)-l]yv = 0. (42) 



Note that we absorb all the difference between GR and MG 
in T] in the above equation, which is the same quantity we 
have shown in the top panels of Fig. [2] Equations ( 32 33 



40 42 I form a set of coupled nonlinear differential equations, 



which govern the evolution of an underdensity in a given 
environment. 



We can now solve the above equations to track the evo- 
lution of a spherical top-hat void. We compare results in our 
chameleon model and in ACDM in Fig. |3j both start from 
the same underdense regions 5sc = —2.76, where 5sc (shell- 
crossing) is the initial density contrast of the void region 
extrapolated to today. This setting of the initial condition 
corresponds to voids that would have just shell-crossed to- 
day in the ACDM universe ( Sheth & van de Weygaert'2004[). 
The mean nonlinear density contrast of those underdense re- 
gions today is (5 = —0.8, so that even without the fifth force 
these are already fairly empty voids. 

The difference between the two models in the void ex- 
pansion history depends on the initial comoving sizes of 
voids R as well as their environment, quantified by (5cnv, the 
initial environment density perturbation linearly extrapo- 
lated to today. Voids in denser environments show a larger 
difference between GR and MG. This is due to the greater 
density contrast realized by an underdensity in a very over- 
dense environment. As seen in Figs, [l] and [2] such contrasts 
in density cause a large change in the scalar field, which in 
turn results in a stronger fifth force. 

In all cases, voids in MG expand faster and grow larger 
than their counterparts in ACDM. The comoving void ra- 
dius would have grown by a factor of 1.7 at shell-crossing in 
GR. However in MG, Fig. [3] shows the same underdensity 
would have grown by a factor of ~ 2 for voids with i? ~ 3 
Mpc/ft in dense environments. The difference between GR 
and MG is at ~ 10% level, and smaller for less dense environ- 
ments. For larger voids the difference becomes smaller, e.g., 
for R ~100 Mpc//i it is at the sub-percent level. Although 
the absolute value of the fifth force is smaller for small voids 
(left panel of Fig. [T]), the gravitational force is correspond- 
ingly smaller due to the decreased integrated mass. As shown 
in Fig. |2] the net effect is that the instantaneous ratio be- 
tween the two is larger (more negative) for smaller voids. 
The void size y^ or r at any given time shows the integrated 
effects of nonzero 77 from all previous times. Thus the radii 
of smaller voids have expanded more beyond their ACDM 
counterparts, which themselves have expanded much more 
than the background. 

While the size of voids shows the cumulative effect of 
gravity, the expansion velocity of each shell responds more 
sensitively to any change of gravity at a given time. The 
bottom panels of Fig. [3] shows the fractional difference of 
the expansion velocity of shells in GR and MG. Indeed, the 
differences in velocity are larger than the differences in sizes. 
For voids of i? ~ 3 Mpc//i, the expansion velocity can be 
10% to 30% faster in MG in over-dense environments. By 
_R ~ 30 Mpc//i, the difference has dropped to a few percent 
in this model. 

Our results suggest that perhaps the best way to look 
for modified gravity is to find voids in overdense environ- 
ments, especially small voids, where we expect the differ- 
ence from GR is maximized. Those voids should be emptier 
due to the relatively strong repulsive fifth force and faster 
expansion of the shells. Moreover, the difference in redshift 
space could be more prominent due to the even larger dif- 
ference in the velocity field. We propose that the clustering 
analysis of tracers of small voids in redshift space could be a 
powerful test of GR. Predictions for this test from A'^-body 
simulations will be presented in a separate paper. 
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R = l Mpc/h 



R=3 Mpc/h 



R = 10 Mpc/h 



R =30 Mpc/h 




Figure 3. Top row: Radius r of an expanding underdensity in units of its initial comoving radius R, as a function of scale factor a. 
Center row: Fractional difference between the radii of such underdensities with identical initial conditions, expanding with and without 
the fifth force. Bottom row: Fractional difference in the velocity. Columns show various values of initial comoving radius, R = f,3,f0 
and 30Mpc//i, from left to right. All panels have an initial underdensity, linearly extrapolated to today, of 5 = —2.76: these are objects 
which in a universe with no fifth force would have just reached the epoch of shell-crossing today. Various values of the exterior density 
are shown, with (5onv decreasing from top to bottom. The largest deviations from GR occur for voids expanding within a larger overdense 
region. 



5 VOID DEFINITION AND STATISTICS 

Having success in following the evolution of a single shell, 
we can now look for a common definition of voids for GR 
and MG. Then we will compare the population of voids in 
both GR and MG statistically by generalizing the excursion 
set approach ( Bond et al.||199T| ). But first we will lay down 
briefly the essential idea of the excursion set theory; more 
details can be found in Appendix [X] 



5.1 Excursion Set Theory 

Assume that the initial local density perturbation filtered 
at a given scale R, 5{x, R) follows a Gaussian distribution, 
and that there is no correlation of 5{x, R) between different 
filter sizes (for correlated 5, see Musso & Sheth (2012|). 
Then we know A.) the distribution can be fully described 
by its variance 5", and B.) when varying the filter size 7? to 
R — dR or equivalently in hierarchical models, S — >■ 5* + dS, 
the increment of 5{x, R) is independent from its previous 
value and should also follow a Gaussian distribution with the 
variance of dS. Thus, S{x, S) is just a Brownian motion with 
'time' variable S. In the spherical collapse model, if a local 
density exceeds a certain barrier 5c, then it will collapse and 



form a virialized halo with all the mass M' enclosed within 
R' by some given time. In the {S, 5)-plane, if we start the 
walk from the origin, walks that cross 5c for the first time at 
S' — a'^{M') correspond to such objects. Walks which cross 
first at smaller values of S form higher mass halos. Therefore, 
the fraction of mass that has collapsed and formed halos 
heavier than M' is the fraction of random walks 5{x,S) that 
have crossed the barrier 5c at S < S' . Alternatively, one can 
calculate the fraction of mass that is incorporated in halos 
at a given range of halo mass [M , M + dM], or equivalently, 
[S, S + dS] at a given redshift z: 



fiS,z)dS- 



D+iO)5c 



exp 



Dl{0)5l 



2Dl{z)S 



dS, 



(43) 



where f{S,z) the first-crossing distribution of the Brown- 
ian motion to the barrier D+{0)5c/D+{z), and D+ is the 
linear growth factor. The first crossing distribution essen- 
tially gives the halo mass function (see Appendix [A|. There 
is equal chance for a random walk to go negative in S. Thus, 
once an appropriate first-crossing barrier for voids, 5^, is 
given, one can also find the void size distribution function 
by the same method. 
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5.2 First crossing barrier for void 

For halos, the first crossing barrier 5c is usually defined as 
the linearly extrapolated initial overdensity at the time of 
collapse, i.e., when the mass shells reach zero radius. This 
time can be calculated using the spherical collapse model. 
Naively one can find the shell-crossing barrier for voids in a 
similar way. The shell at the radius of r of a perturbed spher- 
ical underdense region will expand faster than the shell at 
r' = r 4- Ar, as the enclosed mass within the border shell is 
smaller. Shell-crossing for underdense regions occurs when 
the two shells collide. This occurs at the present day for 
underdense regions with 5bc = —2.76 (the density contrast 
at the initial condition extrapolated to today) for the con- 
cordance ACDM model. Like 5c, 5sc depends on Q.^ and is 
independent of smoothing scale. Moreover, the underdense 
region at shell-crossing happens to be very empty, i.e., its 
nonlinear underdensity is 5 = —0.8. Therefore, 5 = —0.8 
serves nicely as an empirical definition of voids. 

In modified gravity, however, the situation is more com- 
plicated. First, the shell-crossing barrier can depend on the 
environment, simply because the fifth force and hence the 
expansion history of shells depends on the environment. 
Therefore, one may expect voids (likewise halos ( |Li &c Efe^ 
tathiou|2012 l) to form differently depending on the environ- 
ment. Second, even for the same environment, the popula- 
tion of voids may also be different from ACDM, due to the 
size dependence of the force which leads to scale dependence 
of the barrier. 

In chameleon models, the fifth force does speed up the 
expansion of voids (as seen in Fig.[3|, but the shell-crossing 
time usually occurs later than in ACDM with the same ini- 
tial conditions. This is because the effect of the fifth force 
on the relative accelerations of neighboring shells is in the 
opposite direction from gravity. For — 1 < r; < the fifth 
force opposes but does not overcome gravity, so that the 
pull of inner shells on outer ones is reduced, making the crit- 
ical density for shell crossing in chameleon models harder to 
reach. If an observer is riding on the boundary shell, then 
all the nearby shells move closer with time, but more slowly 
than shells feeling only standard gravity. Furthermore, for 
some initial density perturbations, the shell crossing does 
not happen at all. 

Since the epoch of shell crossing can be unreasonably 
late or undefined for these models, it is easier to use em- 
pirical criteria for void formation. We choose 5 = —0.8 as a 
common criteria for the following reasons. A.) it correspond 
to the first-crossing barrier in ACDM, making it easy to 
compare with results from ACDM; B.) Voids with 5 = —0.8 
are indeed very empty, and can be defined by the same way 
in simulations and observations, thus enabling one to make 



direct comparisons. For example, in Hoyle & Vogeley (20041 



and Sutter et al. (20121 they use similar threshold to define 



voids in the 2dFGRS and SDSS galaxy samples. |Pan et aT] 
(20121 also find voids from SDSS7 having similar density 
contrast, 5 < —0.85 at the edges. 

Thus, we use the requirement that the nonlinear density 
constrast today is 5 = -0.8, along with Eqs. ( [32) [33] [io] [42| , 
to solve for the initial underdensity as a function of scale S 
and environment 5env. The resulting void-formation barrier 
is shown in Fig. [4] Unlike ACDM where the crossing barrier 
is fiat, barriers in chameleon models are scale dependent. In 



general, barriers in chameleon models are lower (less nega- 
tive) . Smaller voids have shallower barriers to reach in order 
to form due to the fact that the fifth force in smaller voids is 
relatively stronger (see Fig.[3|, which makes them to expand 
faster. In other words, for reaching the same 5 = —0.8 to- 
day, the necessary initial density contrast for smaller voids is 
smaller (less negative) . The crossing barriers keep decreasing 
(becoming less negative) and steepening with the increase of 
5*. This is very different from the collapsing barrier for halos 
in the same model, where they are leveling off at S ~ 5 ( [lI| 
fc Efstathio"u||2012[ ). This difference is a direct result of the 
fifth force strength upper bound of , which only applies 
to overdensities. 

Fixing void size, the barrier is lower (less negative) 
and steeper for denser environment, where the difference 
from the flat barrier in ACDM is also larger. Therefore the 
difference of void population with ACDM should be more 
prominent in such regions. This environmental dependence 
of crossing barrier is the opposite for halos, where the col- 
lapsing barriers are higher (more positive), and closer to 



the ACDM barrier for denser environment ( Li & Efstathiou 



2012 1. Qualitatively, these two opposite pictures in voids and 



halos can be understood by the same reasoning, i.e., for voids 
or halos of the same mass given pin (the mean density in the 
void or halo region), the strength of the fifth force is larger 
for larger differences between pin and the background den- 
sity outside the perturbed region, pout - For voids *C pout 5 
er pout means pout — Pin! is larger and hence a larger 
fifth force, while for halos pin > Pout, a larger pout means 
I Pout — Pin I is smaller therefore a smaller fifth force. 

5.3 Moving environment approximation 

In calculating the void barriers in the previous section, the 
environment was specified only by its linear density pertur- 
bation, (5cnv. In order to derive the first-crossing distribu- 
tions and other void statistics it is necessary also to specify 
an environment length scale. 

In treating spherical collapse in chameleon models, [LI] 
I&: Efstathiou| ( p012[ | used an environment scale of 8 Mpc//i 
for halos of every size. Such a fixed-environment scale works 
well for halos, since the range of interesting virial radii is 
fairly small, ~ 0.1 — 1 Mpc//i. Furthermore, since through- 
out collapse the proto-halo is always shrinking, there is little 
worry of its size becoming comparable to the environment 
scale. On the other hand, the interesting void sizes we are 
considering range from ~ 1 — 30 Mpc//i, and each will ex- 
pand beyond its initial comoving radius by a factor 1.7 at 
formation. Thus we need to consider the definition of the 
environment more carefully. 

First, the scale of the environment should at least be 
larger than the final size of the void. Secondly, it should also 
be large enough so that the scalar field in the environment 
has space to settle to its minimum. This is to guarantee that 
the boundary condition Eq. ( |29[ ) for the scalar field profile 
equation holds. Third, it cannot be too large because this 
would simply mean using a value very close to the cosmo- 
logical density p,„ for all void environments. Bearing these 
considerations in mind, we introduce a moving environment 
approximation, in which the initial environment scale is a 
function of the initial void scale, specifically -Rcnv = 5 R. 

We notice that in the moving environment approxima- 
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tion, the expanding void shell and collapsing environment 
shell may cross for voids in very overdense environments. 
Therefore we also calculate the first-crossing distributions 
with a large fixed-environment scale of i?cnv ~ 75Mpc//i, 
so that the environment shell begins its collapse much far- 
ther from the void shell. The difference between the two 
approximations is less than 10% for the void scales of obser- 
vational interest, i.e., those ~ lMpc//i and larger; details 
of the comparison can be found in Appendix [B] It follows 
that the results for choices of -Rcnv > 5 R are also less than 
10%, since such environment scales are between our fiducial 
choice Ranv = 5 7? and the fixed environment scale. This level 
of difference, as we will see later, is negligible compared to 
the difference between GR and MG that we are considering. 
Thus, our main conclusions are insensitive to the definition 
of environment. 



5.4 Conditional first-crossing distributions 

5.4-1 Unconditional First Crossing of a Moving Barrier 

The distribution of the first crossing of a general barrier 
by a Brownian motion has no analytic solutions except for 
some simple barriers, e.g., flat ( Bond et al.|199r I and linear 
( |Sheth|1998 Sheth fc Tormen|2002| ). Unfortunately neither 
of these is a good approximation to our barriers in Fig. [4] 
As a result, we follow Zhang & Hui (20061 and numerically 
compute this distribution. We briefly review their method 
for completeness. 



Denote the unconditional probability that a Brownian 
motion starting off at zero hits the barrier b{S) > for 
the first time in [S,S + dS] by f{S)dS. Then, /(S), the 
probability density, satisfies the following integral equation 



fiS) 



giS)+ / dS'f{S')h{S,S'), 



(44) 



in which 
9iS) 
h{S, S') 



'dS 



^db _ b-b' 



dS S- 



P{b, S), 

P{b-b',S - S'), (45) 



where for brevity we have suppressed the 5'-dependence of 
b{S) and used b' = b{S') and 

2S 



P{S, S)dS 



: exp 



dS. 



(46) 



Equation ( 44 1 can be solved numerically on an equally- 
spaced mesh in S: Si = iAS with i = 0,1,--- ,N and 
AS = S/N. The solution is ( [Zhang fc Huil[20M l 



/o 
/i 



50 = U, 
(1-Ai,i)-V, 



(47) 



(1- Ai, 



5, + ^/,(A i.j + Aij- + l) 



where we have used /; = f{Si) and similarly for gi to lighten 
the notation, and defined 



A, ^ 



AS'. 



h I Si, Sj 



AS 
2 



(48) 



We have checked that our numerical solution matches the 
analytic solution for the fiat-barrier crossing problem. 



5.4-2 Conditional First Crossing of a Moving Barrier 

The unconditional first crossing distribution, which relates 
directly to the void size distribution function in the ACDM 
model, is not particularly interesting in the chameleon 
model. This is because spherical underdensities in differ- 
ent environments will follow different evolution paths. If 
it is in the environment specified by (Scnv, 5cnv), then 
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(5'onv,5env) should be the starting point of the Brown- 
ian motion trajectory. In other words, we actually require 
the distribution conditional on the trajectory passing Scm 

at S —— Scnv 5 

we write this first-crossing distribution as 
f{S,dviS,Scnv) I Scnv,(5cnv), showing exphcitly the Scnv de- 
pendence of 5v The numerical algorithm to calculate the 
conditional first crossing probability is a simple generaliza- 
tion of the one used above to compute the unconditional 



first crossing probability (Parfrey et al. 20111 and is not 
presented in detail here. 

Note that the preceding algorithm assumes the barrier 
b{S) > 0, while our void-formation barriers are strictly neg- 
ative. However, if solving the problem by a Monte Carlo 
method we could note that the resulting first-crossing dis- 
tribution is invariant under reflecting the Gaussian ran- 
dom walks about S — (since each step of each walk is 
equally likely to move to higher or lower S). Thus, we can 
solve the distributions for our negative barriers by using 
b{S) = |5v(S')| in the above algorithm. 

Furthermore, the preceding algorithm describes the 
calculation of the first-crossing probability for the fixed- 
environment approximation, in which a single starting point 
(iSenv, <5env) for a giveu barrier 6viS,Senv) is sufficient. To 
implement the moving environment approximation we cal- 
culate a new first-crossing distribution for each underden- 
sity scale S, where the walk starts at Scnv{Rcnv) and 
Renv = 5R{S) as described in Sec. 5.3 Our final re- 
sult for the conditional first-crossing probability is then 
f{S,5v{S,5cnv) I Scnv{S),5cnv), wherc the dependence of 
Scnv on S is written explicitly. 

In the special case where the barrier is flat, Sv{S, Scnv) ~ 
Ssc, f{S,5v{S,5cnv) I Scnv{S),Scnv) IS kuowu analytically as 



/ 



I ^sc ^cnv I 



27r(S'- 



^3/2 



exp 



2{S- Sc, 



(49) 



where again Scnv = 5oiiv(S'), so that in the next section we 
compare first-crossing distributions for GR and MG both 
calculated using the same moving environment scale. 



5.4.3 Results 

Figure [5] shows the first-crossing distribution of voids in 
different environments. In general, we find all voids today 
with radii -R/ > 1 Mpc//i are more numerous in chameleon 
models, for all environments. This difference from ACDM is 
larger for overdense environments. This is a consequence of 
previous results of this paper, namely that the fifth force is 
relatively stronger for denser environments. 

Next, consider fixing the environment density. In this 
case, the fractional difference of the number density between 
chameleon models and GR tends to be greater for larger 
voids (larger u = S^c/S or smaller S), as indicated by the 
increase of A/// with v in the figure. For example, in the 
environment of Scnv = 0.8, voids with Rf = 5 Mpc/h may 
be 2 to 3 times more common than those in ACDM, and 
10 times more for Rf = 25 Mpc/h. This difference may 
seem surprisingly large, but such a case may be too rare to 
be observed. If one smooths the initial density field with a 
filter size much greater than 7i = 15 Mpc/h (correspond- 
ing to Rf —25 Mpc/h), the probability distribution of the 
overdensity will be a narrow Gaussian with zero mean. The 



chance of having a linearly-extrapolated Scnv ~ 0.8 should 
be very low; the odds of such an environment developing 
voids of Rf = 25 Mpc/h or larger with S — —0.8 will be 
even less. Therefore, it might be difficult to find large voids 
in very overdense environments, where the predicted differ- 
ence between models is expected to be larger. In reality, most 
large-scale environments are very close to the cosmic mean, 
i.e., Scnv ~ 0. In this case, the difference between models 
indicated by A/// is less extreme but still very significant, 
being ^ 100% for Rf = 5 Mpc/h and ^ 300% for Rf = 25 
Mpc/h. We shall see in the next subsection that this differ- 
ence is indeed close to the case where the average over all 
environments is taken. 

The environmental dependence of model differences in 
the conditional first crossing distribution of voids is just the 
opposite as that for halos for reasons we have explained in 
Sec. 5.2 The halo mass function (Li & Efstathiou 2012) is 



found to differ more from its ACDM counterpart in under- 
dense environments. 

The fact that A/// is larger for larger voids might 
seem counter-intuitive, as we have shown that the relative 
strength of the fifth force is smaller for larger voids, hence 
the difference in their expansion velocities and sizes today 
are relatively smaller. However, the difference in the number 
density of voids is also related to the shape of the void size 
distribution function. Consider that / is a very steep func- 
tion of when 1/ is large. A small increment in Rf or v can 
therefore lead to a relatively large change in /. 

In principle, if /mg is larger than /gr for large voids, 
the opposite should be true for small voids, namely the 
abundance of small voids will be lower in chameleon models. 
This is expected from the normalization of the first-crossing 
probability. Picturing this in the excursion set theory, in 
chameleon theories Brownian motions are likely to cross the 
barrier at a slightly earlier 'time', i.e. small 5", corresponding 
to large voids. Correspondingly, the probability of a Brow- 
nian motion to survive for longer and cross the barrier at 
large S is reduced - voids of smaller sizes are (relatively) 
rarer than in ACDM. Therefore, the solid and dashed lines 
in Fig. [5] will cross each other at some small u that is not 
plotted, namely the abundance of small voids can be lower 
in chameleon models. In fact, such a crossing point is also 
expected for halos, which has been shown to be at S < 10 



for the environments under consideration (Li & Efstathiou 
|2012[ ). For voids, the crossing points are found to appear at 
much larger S. This is likely due to the halo barriers level- 
ing off at 5 ~ 5, while the void barriers continue to steepen 
towards larger S. 

In real observations, one needs to have tracers like 
galaxies or galaxy clusters to define void walls. If the size 
of the void is comparable to that of the tracers, then the 
walls will be lumpy. Voids with radii comparable or smaller 
than the typical size of virialized objects are therefore not 
well defined and of little interest. We do not show results 
deeply into this regime. In the range of empirical interest, 
we only see the lines of /mg and /gr crossing each other 
for the case of Scnv = —2.4 at ii ~ 1 Mpc/fe, which should 
be a rare situation. Thus, for denser environments we al- 
ways expect to find more voids in chameleon models at all 
empirically meaningful sizes. 

The environmental dependence of the differences be- 
tween models may provide useful guidelines for testing grav- 
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Figure 5. Left: First crossing distribution functions for different environments as indicated by (5cnv Solid lines are in chameleon 
cosmology with our fiducial model parameters. The top x-axis labels the corresponding final void radius when the density contrast of a 
void reaches 5 = —0.8. In bottom x-axis, Ssc = —2.76 is the shell-crossing barrier for voids in ACDM, and S = a'^(M) is the variance 
of a spherical top-hat region. Right: Fractional differences of the first crossing distributions between GR and chameleon cosmology for 
different environments. 



ity. In overdense environments, one may want to look at the 
statistics of large voids as the difference with ACDM may 
be larger, while in underdense regions, the difference in halo 
population may be larger therefore halo number densities 
may be more interesting to analyze. We summarize these 
two cases as void-in-cloud and cloud-in-void. However, both 
of these two cases are relatively uncommon to find in the 
real universe so that the statistics may be poor. In this case, 
using most of the observed volume could provide better con- 
straints since the sample of voids and halos would be larger. 
It is therefore interesting to determine the overall difference 
between models once we average over all different environ- 
ments. 



5.5 Environment-averaged first-crossing 

To get the average first crossing distribution of the mov- 
ing barrier, we must integrate over all environments. The 
distribution of 5onv, denoted as g(5env, ^c, S'onv), in which 
5c is the critical overdensity for spherical collapse in the 
ACDM model, is simply the probability that the Brownian 
motion passes (5cnv at Scnv and never exceeds 5^ for 5* < Scm 
(because otherwise the environment itself has collapsed al- 
ready). This has been derived by Bond et al. (19911: 



Q(5cnv, 5c, iS'cnv) 



: exp 



: exp 



(5c 



25c 



2S, 



,(50) 



for 5env ^ 5c and otherwise. Again, we have Scnv = 
Senv{S), so that the distribution q changes for each void 
size. For smaller smoothing length (larger Scnv), the pdf of 



5cnv is wider so that the very overdense and very underdense 
environments are more likely to be sampled. 

Then the environment-averaged first crossing distribu- 
tion will be 

/avg(S) = f qX f{S, 5^{S, 5cnv) ! 5cnv(5'), 5cnv) d^M) 

The environment-averaged first-crossing distribution and 
void volume function are related by 



dn 
dV 



dV 



M 



/avg(S') 



ds; 

dV 



dV, 



where V is the final volume of the void given by 

y = — X {i.nf. 

pm 



(52) 



(53) 



The factor of 1.71 results from our void formation criteria 
of nonlinear density 5 = —0.8. By mass conservation, such 
an underdensity which was originally at the cosmic mean 
has grown to 1.71 times its initial comoving radius. 

The left and right panels of Fig. [6] show the 
environment-averaged first-crossing distribution and the 
corresponding void volume function, respectively. Compar- 
ing the environment-averaged void distribution functions be- 
tween our fiducial chameleon model and ACDM, we find the 
fractional difference in the number density of voids between 
the two models increases with void size. At Rf ~ 25, one 
may expect to find 2 to 3 times more voids in chameleon 
models, and such a difference will keep increasing for larger 
voids. This level of difference in the void population is much 
greater than that in halos, where the difference of mass func- 



tion is found to be no more than 20% ( Li & Efstathiou 20121 



a factor of 10 times smaller difference. The boost of probabil- 
ity for having large voids in chameleon models has interest- 




Figure 6. Left panel: Compares the averaged first crossing distribution functions between chameleon model (solid line) and GR (dashed 
line). The fractional difference is shown in the bottom panel. Right panel: Void volume distribution functions and their fractional 
difference. The difference in the number density of voids between the two models increases monotonically with void size. 



ing implications for observation, thus serving as a powerful 
test of gravity theories. Given a finite survey volume, one 
can simply count the number of voids greater than a certain 
radius, e.g., i?/ > 25 Mpc/ft to find out the number density 
of them and then compare it with different models. 



changes due to the steepness of the volume function there: 
for a = 10~^ the deviation of the volume function from the 
GR result grows by a factor of 10 in moving from 27^ = 1/3 
to 27^ = 1. 



5.6 Theory Variations 

Up to this point, we have only shown results for our fiducial 
chameleon theory, with parameters a — W^'' and 27^ — 1/3. 
Figure [T] shows the effect of varying these two parameters 
on the volume function, dn/dV. Focusing on the leftmost 
panels we see the models with 27^ = 1/3, which correspond 
most closely to the f{R) class of theories. In moving from 
a = 10~^ to 10~^ the fractional difference changes by a 
factor ~ 3 for smaU voids (V = 7 x 10^ (Mpc/hf) and by 
~ 25 for voids two orders of magnitude larger (V = 7 x 
(Mpc/hf). 

A direct comparison of this chameleon theory with 
f{R) models is not possible, but comparing the compton 
wavelengths can give some idea of the differences. For the 



f{R) model of Hu & Sawicki (20071, the compton wave- 
length in the background density today is ~ 3 Mpc/h for 
I /ho I ~ 10~®. Our fiducial model has a longer compton 
wavelength: for 27^ = 1/3 we have A ~ 2VlO«a. Thus for 
a = 10"^ A ~ 20 Mpc/h. 

As it is also interesting to put constraints on the cou- 



pling 27^ (e.g., (Jain et al. 2012 1), we show such variations 



in the center and right panels of Fig. [7] As we expect, for 
fixed a the deviations are much larger for stronger couplings. 
Again the largest, rarest voids are most sensitive to these 



The comparison to results for the excursion-set mass 
function highlights the promise of using voids to constrain 
modified gravity. Consider the case in Fig.[7|with the small- 
est deviations from GR, a = 10^' and 27^ = 1/3. The 
fractional difference in the volume function is 30-60% over 
at least two decades in void volume. The deviation of the 
mass function predicted by this model peaks at 5% for halo 
masses ~ W^^Mq/H, falling quickly for smaller and larger 



halos ( Li & Efstathiou 2012 1. Thus, if the difference between 



models is integrated over the entire range of halo and void 
number densities, the total constraining power of the void 
statistics will be much greater. This larger difference in the 
void statistics is a result of several effects: A.) the upper 
bound in the ratio of gravity and the fifth force does not ap- 
ply to underdensities and B.) the crossing point of GR and 
MG first-crossing distributions expected due to the normal- 
ization of the distribution occurs for voids which are too 
small to be empirically relevant. Thus, the MG void volume 
function shows large deviations at all void sizes. 



In Appendix [C] we discuss the effect of varying a and 
7 on the conditional first-crossing distributions, i.e., before 
the environment averaging is carried out. 




Figure 7. Left panels: The void differential volume function for our chameleon model with 27^ = 1/3 (solid lines) compared to GR 
(dashed). Various values of a are shown, ranging over 10~^, lO"'^, and 10~^, from top to bottom. The lower panel shows the fractional 
difference from the GR result. Center panels: The same, but for coupling 27^ = 1/2. Right panels: The same, but for coupling 27^ = 1. 
Even for a = 10~^, 27^ = 1/3 where the deviation is weakest, it is above 30% for all empirically interesting void sizes. 



6 CONCLUSION AND DISCUSSION 

We have explored the physics of the fifth force in voids 
for chameleon models and applied it to understand the im- 
pact on void properties. In scalar-tensor theories, such as 
chameleon MG, the smooth part of the scalar field is the 
source of the cosmological constant, known to act like a re- 
pulsive force. This is common in both a ACDM universe 
and a chameleon universe. The coupling of the scalar field 
to mass density causes an additional spatial fluctuation of 
gravity, i.e., the fifth force. This is the only difference for 
void evolution between chameleon and ACDM models. The 
evolution of voids in MG is affected by the Newtonian force, 
the dark energy force and the fifth force. 

The following interesting features are found in com- 
parison to a ACDM universe, some of which may be used 
to test gravity in laboratory experiments and observational 
data, or to guide more precise predictions from cosmological 
Ai'-body simulations. 

1. ) The fifth force in voids is a type of 'anti-gravity'. 
It points outwards from the center of the void, opposite to 
the direction of normal gravity. This is because the slope of 
the scalar field profile is negative in voids. 

2. ) In principle, the amplitude of the fifth force can 
be very large in voids. In halos, its magnitude can be no 
more than 27^ (1/3 in our fiducial model) of normal gravity. 
Due to the breakdown of Birkhoff's theorem, the scalar 
field profile and hence the fifth force are functions of the 



matter density inside and outside the void region, as well 
as its size. 

This property leads us to suggest a possible laboratory 
test of gravity using a vacuum chamber. To create a 
chamber inside of which the fifth force is dominant, it 
should have a thick chamber wall made of high density 
material. This is to enlarge the density contrast between the 
chamber interior and the wall so that fifth force strength 
is maximized. The wall needs to be thick to have enough 
space for the scalar field to reach its minimum in the wall. 
Walls of the chamber and test particles (detectors) in the 
chamber should feel the flfth force pushing outwards, but 
very little gravity (as long as the chamber is close to real 
vacuum) . In the neighborhood of the earth, the background 
density is non-zero. There is dark matter from the Milky 
Way halo, and maybe some baryonic dust; these two 
should contribute a haze of mass density inside any vacuum 
chamber. This may set the limit for the density contrast and 
the amplitude of the flfth force. Furthermore, although the 
ratio -FsZ-Fn may be large in this case, we know that _Fn is 
quite small in the chamber, so that the large ratio does not 
necessarily imply a large flfth force. The effect of the dark 
energy force also needs to be accounted for. We leave the 
quantitative investigation of this experiment to future work. 

3.) Driven by the additional fifth force, individual 
voids expand faster and grow larger than their ACDM 
counterparts. The fractional difference in void radius and 
expansion velocity is larger for small voids in overdense 
environment (void-in-cloud) , at the level of < 10% and 
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20 — 30% respectively, for voids of a few Mpc/fe. For the 
same reason, voids of the same size should be emptier in 
chameleon models. This leads to interesting observational 
consequences. A.) In redshift space, due to the faster 
expansion of voids, a small void-in-cloud may be more 
elongated along the line-of-sight due to redshift space 
distortion. B.) Void profiles may be steeper as voids empty 
themselves more quickly, as has been shown in |Martino fc| 
Sheth (20091. We plan to investigate both of the above by 
stacking voids in simulations. Recent work has shown that 
the lensing signal from stacked voids in future surveys will 



Local Void and piling up matter on its outskirts" (Peebles 
& Nusser 20101. Voids in chameleon models seem to 



coincide qualitatively with these observations. However, the 
complexity of galaxy formation, especially its dependence 
on environment, is a hard barrier to overcome before any 
conclusive results can be drawn. 

We note that our results for void statistics should be 
qualitatively similar in other models with chameleon 
screening, such as ,f{B) ( Hu fc Sawicki|2007 1. Furthermore, 



20121. This may provide a complimentary probe to void 



while symmetron ( Hinterbichler & Khoury 2010 Hinter- 



provide information on their radial profile (Krause et al. bichler et al.||2011 l and environmentally-dependent dilaton 



statistics for distinguishing between gravity models. 

4.) For individual voids, the largest difference between 
GR and MG is found in void-in-cloud systems, while for 
voids statistics, the large voids differ more. The fractional 
difference in the number density of voids increases with size 
and is --^10 times larger than the corresponding difference 
for halos. The chance of having voids with 5 ~ —0.8 with 
7? ~ 25 Mpc/h is 2.5 times larger than in ACDM. A 
conceptually simple observational test would be to count 
the number of very large halos in a volume limited sample, 
and find out the probability for that count to occur within 
different gravity models. 

In fact, the detection of the CMB Cold Spot in the 
WMAP data, if interpreted as the ISW signal, has already 
imposed a constraint on this probability. The size of the void 
in the large-scale structure needed to generate the size and 
amplitude of the Cold Spot is estimated to be at the or- 
der of 100 Mpc//i in radius, which may not be consistent 
with a ACDM universe (e.g . |Cruz et al. |2005| 



Inoue & Silk 

2006[ [Rudnick et al.|[2007[ [Masina fc Notari|r2009| |Cai et 
al.|2010[ |. Similarly, the detected ISW signal from the stack- 



ing of 4-deg^-size regions of the CMB corresponding to the 
SDSS super clusters and super voids is found to be 2 — So- 



higher than that expected in a ACDM universe ( Granett et 
al.|2008 Papai et al.'2011) . Recent work has shown that the 



abundance of the largest voids in ACDM simulations may be 
too small to match observations ( Tavasoli et al.|["2012 |. All 
of these discrepancies, if confirmed, seem to indicate that 
very large structures in the universe are perhaps larger and 
more abundant than expected in a ACDM universe. The fact 
that the abundance of large voids in modify gravity can be 
much greater than in ACDM suggest that modify gravity 
can somewhat release the tension imposed by those obser- 
vations, but precise quantitative predictions are beyond the 
reach of the spherical collapse model and excursion set the- 
ory. 

Intriguingly, there are also observations suggesting 
that galaxies are less common in low density regions than 
expected in the standard cosmology (e.g. /Tikhonov & 



Klypin|[2009| ). The Local Void (within the radius of 1-8 
Mpc from the center of the local group) also seems far too 



empty based on the galaxy number density (e.g. Peebles 
20011 iTikhonov fc Klypin"2009| iPeebles fc Nusser||2010l), 



but see [Tinker fc Conroy (2009) for a different view. There 
is also an unexpected presence of large galaxies on the 
outskirts of the Local Void ( Peebles fc Nusser|2010 1 . "These 
problems would be eased if structure grew more rapidly 
than in the standard theory, more completely emptying the 



(Damour fc Polyakov 1994 Brax et aLjpOlO I theories rely 
on conceptually different mechanisms to screen the fifth 
force, the qualitative picture of Fig. [l] is unchanged. The 
minimum of the symmetron and dilaton fields will again be 
higher inside an underdensity than outside, thus leading 
to a repulsive fifth force which will aid the dark energy in 
speeding up void growth. 

Caveats: 

Throughout the paper, we employ the spherical collapse 
model and excursion set theory for studying the evolution of 
individual voids and their distribution functions. However: 

1. Voids in the real universe are not perfectly spherical, 
(e.g. Shandarin et al.|2006 |. 

2. The excursion set theory for voids may not be able to 
match precisely voids found from simulation or observation. 
There are obvious reasons for this:. A.) It has been noticed 
that the total volume of voids given by this model exceeds 



that of the universe ( Sheth fc van de Weygaert||2004i ). This 
is certainly not physical. One obvious reason is that some 
'voids' may be embedded in overdense regions whose den- 
sity reaches the collapsing barrier. This is the void-in-cloud 
problem, which is more acute for small voids. Accounting 
for it can resolved the problem to some extent, but not fully 



( Sheth fc van de Weygaert|2004[|Paranjape et al.|2012[ ). An- 
other reason is that there is an underlying assumption that 
voids can expand forever, which is also unphysical. The ex- 
panding walls of voids will certainly meet their neighbors 
and cross each other. This is probably more complicated to 
fix and is beyond the scope of this paper. 

3. Our results are for voids in the dark matter distribu- 
tion, whereas observed voids must be defined with respect 
to galaxies. The excursion set theory of the void population 
has been extended to these more empirically relevant voids 



by Furlanetto & Piran ( 2006 1 



In this paper, we are mostly comparing the difference 
between two models rather than the accuracy of each model 
itself. Thus, these well-known limitations of the basic ex- 
cursion set theory of voids may affect MG and ACDM in 
roughly the same way, leaving the difference mostly unaf- 
fected. We therefore neglect these problems, and leave the 
calibration of the theory to simulation for future work. 
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APPENDIX A: EXCURSION SET THEORY 

It is widely accepted that the large-scale structure (LSS) 
in the Universe has developed hierarchically through grav- 
itational instability. The excursion sets (regions where the 
matter density exceeds some threshold when filtered on a 
suitable scale) generally correspond to sites of formation of 
virialised structures (jSchaeffer & Silk 1988| Cole & Kaiser 
19881 [Co ie & Kaisers 1989; E fstathiou et al.||1988| |Efs- 
tathiou fc Rees 1988, ,Narayan ^ White||1987| fcarlbergT 
Couchman|1988] 



The filtered, or smoothed, matter density perturbation 
field (5(x, R), is given by 



5(x,i?) 



VK(|x-yi;i?)5(y)dV, 



VK(fc;E)<5ke*"d^k, 



(Al) 



where W{r; R) is a filter, or window function, with radius 
R, and W{k;R) its Fourier transform; (5(x) = p(x)/p— 1 is 
the true, unsmoothed, density perturbation field and 5k its 
Fourier transform; we will always use an overbar to denote 
background quantities. 

As usual, we assume that the initial density pertur- 
bation field (5(x) is Gaussian and specified by its power 
spectrum P{k). The root-mean-squared (rms) fiuctuation of 
mass in the smoothing window is given by 



S{R) 



HR) = {S\^;R)) 



P{k)W{k;R)<rk. (A2) 



Note that, given the power spectrum P{k), S, R and M are 
equivalent measures of the scale of a spherical perturbation 
and they will be used interchangeablly below. 

If W{k; R) is chosen to be a sharp filter in k-space, then 
the increment of S{x; R) as R ^ R — 6R or equivalently 
5* — >■ S + SS comes from only the extra higher-fc modes 



of the density perturbation (see Eq. (All). The absence of 



correlation between these different wavenumbers means that 
the increment of (5(x; R) is independent of its previous value. 
It is also a Gaussian field, with zero mean and variance SS. 
Thus, considering 5* as a 'time' variable, we find that 5(x; S) 
can be described by a Brownian motion. 

The probability distribution of 5(x; R) is a Gaussian 



P{5, S)dS 



: exp 



5^ 
'2S 



dS. 



(A3) 



In an Einstein-de Sitter or a ACDM universe, the linear 
growth of initial density perturbations is scale-independent, 
so that 5(x) and cr{R) — \/S grow in the same manner, and 
as a result the density field will remain Gaussian while it is 
linear. Following the standard literature, hereafter we shall 
use (5(x; R) to denote the initial smoothed density perturba- 
tion extrapolated to the present time using linear perturba- 
tion theory, and the same for a ot S. 

In the standard cold dark matter scenario, the initial 
smoothed densities which, extrapolated to the present time, 
equal (exceed) Sc correspond to regions where virialised dark 
matter halos have formed today (earlier). In an Einstein-de 



Sitter universe Sc is a constant, while in a ACDM universe it 
depends on the matter density fim. In neither case does Sc 
depend on the size of (or equivalently the mass enclosed in) 
the smoothed overdensity, or the environment surrounding 
the overdensity. 

As a result, to see if a spherical region with initial radius 
R has collapsed to virialised objects today or lives in some 
larger region which has collapsed earlier, we only need to see 
whether 5(x; ^ R) ^ Sc- Put another way, the fraction of 
the total mass that is incorporated in virialised dark matter 
halos heavier than M = ^nR^pi is just the fraction of the 
Brownian motion trajectories 5(x; S) which have crossed the 
constant barrier Sc by the 'time' S = S{R), which is given 
by [Bond et"ar]|T99Ti 



F{M,z) = 



D+(0) 



d5. 



(A4) 



where the lower limit of the integral is 



P+(0) 
D+(z) 



Sc, because if 



a virialised object formed at redshift z, then its correspond- 
ing initial smoothed density linearly extrapolated to z is Sc, 



while extrapolated to today it is 



Sc with -D+(z) being 



the linear growth factor at 2. In Einstein-de Sitter cosmology 
D+{z) oc (1 + z)~^ and this quantity becomes (1 -|- z)Sc. 

Alternatively, one can say that the fraction of the total 
mass that is incorporated in halos, the radii of which fall in 
[R, R + SR] (or equally [S, S + SS]) and which collapse at 
z = Zf is given by 



f{S,Zf)dS^ 



1 D+{0)Sc 
/2^D+{zf)S 



exp 



2Dl{zf)S 



dS, (A5) 



where f{S) the distribution of the first-crossing time of the 
Brownian motion to the barrier D+{z = 0)Sc/ D+[z = Zf). 
Once this is obtained, one can compute the halo mass func- 
tion observed at z/ as 



'^"(^)dM = ^4^/(5)dS. 



dM 



M 



(A6) 



Other observables, such as the dark matter halo bias (Mo 
fc White! 19 9"6| or merger history (Lacey & Cole 1993), can 
be computed with certain straightforward generalizations of 
the theory. 



APPENDIX B: COMPARING FIXED- AND 
MOVING-ENVIRONMENT MODELS 

In the main text we use a moving-environment approxima- 
tion in which the smoothing scale of the environment is a 
function of the void scale, specifically i?cnv ~ 5R. However 
we have also checked the effect of using a fixed-environment 
approximation to calculate the fifth force. We compare the 
effect of the approximations on the environment-averaged 
first-crossing distribution in Fig. |Bl[ for a fixed-environment 
scale of -Renv = 75Mpc//i, corresponding to Sonv ~ 0.01. 
The differences are below 10% for v ^ 1, corresponding to 
final void radii Rf > lMpc//i. Thus, throughout the range 
of observable void sizes our conclusions are fairly insensitive 
to the precise approximation used to calculate the environ- 
mental effect of the fifth force. 
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Figure Bl. Upper panel: Environment-averaged first-crossing distribution of voids with (solid lines) and without (dashed lines) the 
fifth force. The higher solid and higher dashed lines show results for the fixed-environment approximation, while the lower pair show the 
moving environment approximation. Lower panel: Fractional difference of fixed- from moving-environment approximation for modified 
gravity (solid) and GR (dashed). For f ~ 1 and larger, the observational range of interest, the difference is below 10%. 



APPENDIX C: THEORY VARIATIONS 

Figure [CT] shows the results for the conditional first-crossing 
distributions, for various parameter values. The results for 
any individual panel are qualitatively very similar to those 
for our fiducial model, a = 10"'^,7 = 1/3. The mam excep- 
tion is for the a = 10~^ theories in very underdense environ- 
ments, 5onv ~ —2.4. Here the random walk begins close to a 
barrier which is itself very near to the ACDM barrier. This 
situation shows that the monotonic increase of the deviation 
with void size is not universal. 

In general, larger values of a allow for much greater 
variation in the conditional first-crossing distributions for 
various environments. Variations in the coupling 27^ cause 
less variation between the different environments. Finally, 
although it is more clearly seen after the environment aver- 
aging (Fig. [7]), larger variations of the distribution with 27^^ 
occur for larger a values. 
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